BONN-TH-2007-05 



Focus Point SUSY at the LHC Revisited 



Siba Prasad Das 1 

Physikalisches Institut der Universitdt Bonn, Nufiallee 12, D-53115 Bonn, Germany 

Amitava Datta 2 

Department of Physics, Jadavpur University, Kolkata 700 032, India 
CN ! Monoranjan Guchait 3 

Tata Institute of Fundamental Research, Homi Bhaba Road, Mumbai 400 005, India 
Manas Maity 4 and Siddhartha Mukherjee 5 
Department of Physics, Visva-Bharati, Santiniketan 731 325, India 



Abstract 



< 
<N 

^p" 1 . The estimation of the backgrounds for gluino signals in focus point supersym- 

metry is extended by including the backgrounds from the production of four third 

generation quarks in the analysis. We find that these backgrounds are negligible 

if one uses the strong selection criteria proposed in the literature (including this 

, analysis) for heavy gluino searches. Softer selection criteria often recommended for 

lighter gluino searches yield backgrounds which are small but numerically signifi- 
ed} ' 

cant. We have also repeated the more conventional background calculations and 
compared our results with the other groups. We find that the size of the total resid- 
ual background estimated by different groups using different event generators and 
hard kinematical cuts agree approximately. In view of the theoretical uncertainties 
in the leading order signal and background cross sections mainly due to the choice of 
the QCD scale, the gluino mass reach at the LHC cannot be pinpointed. However, 
requiring a signal with > 3 tagged b-jets ( instead of the standard choice of > 2 ) it 
is shown that gluino masses close to 2 TeV can be probed at the LHC for a range 
of reasonable choices of the QCD scale for an integrated luminosity of 300 fb _1 . 
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1 Introduction 



Supersymmetry (SUSY) is one of the most elegant extensions of the Standard Model 
(SM) which naturally stabilizes the Higgs boson mass even in the presence of a very en- 
ergy scale like the grand unification scale (M G ) [1]. This occurs due to the cancellation of 
the quadratic divergences among loop diagrams involving SM particles and their super- 
partners (sparticles). This cancellation, however, looses much of its attractive naturalness 
and smacks of fine tuning if the sparticles are much heavier than the corresponding par- 
ticles. Thus SUSY models with sparticle masses much larger than 1 TeV are usually 
disfavored. 

On the other hand it is well known that if the scalar superpartners of the fermions 
belonging to the first two generations in the SM are very heavy with masses in the multi- 
TeV region then several problems can be avoided without further unnatural adjustment 
of the parameters. Firstly in the most general minimal supersymmetric extension of the 
Standard Model (MSSM) there are unacceptably large loop induced flavor changing neu- 
tral currents unless the masses of the first two generations of squark are assumed to be 
degenerate to a very good approximation without a symmetry underlying this degener- 
acy. This SUSY flavor problem becomes less severe if the above loops are suppressed 
by the masses of the heavy scalars. Secondly SUSY models with complex soft breaking 
parameters have potentially dangerous large loop induced CP violating contributions to 
the electric dipole moments of various particles in general. This SUSY CP problem is also 
tamed if the scalars are very heavy. Therefore models with very heavy scalars without 
sacrificing the naturalness criterion are indeed welcome. 

All the above attractive features are realized in the focus point (FP)/ hyperbolic 
branch region [2, 3] of the minimal supergravity(mSUGRA) model[4]. The SUSY flavor 
and CP problems become less severe [5]. The magnitude of /a, the higgsino mass parame- 
ter, which is fixed by the radiative electroweak symmetry breaking ( REWSB ) condition, 
also turns out to be ~ mz as required by naturalness [2, 6]. 

The low value of ii thus obtained has another important implication. The lightest 
supersymmetric particle (LSP) which is the lightest neutralino (xi) over most of the 
parameter space, is usually a mixture of the electroweak gauginos and Higgsinos. However, 
for small ji the higgsino component in Xi dominates. In this case Xi veT J efficiently 
annihilates into gauge boson pairs leaving open the attractive possibility that LSP could 
be a viable dark matter candidate [7, 8, 9, 10]. In fact it is now widely acknowledged that 
the focus point region of the mSUGRA model leads to a dark matter density compatible 
with the latest WMAP data [11]. Thanks to the above interesting features the collider 
signatures of focus point supersymmetry have attracted wide attention[12, 13, 14, 15, 16, 
17, 18, 19]. 

Although the scalar quarks are rather heavy in this model, naturalness arguments 
imply that the gluinos must be necessarily light [2]. Thus the gluinos are the lightest 
strongly interacting sparticles within the striking range of the LHC. Inspite of the fact 
that the squarks belonging to the first two generations are indeed very heavy, the lighter 
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squark mass eigenstates belonging to the third generation, though significantly heavier 
than the gluino, may be relatively light at the weak scale M w . This is partly due to the 
renormalization group (RG) evolution from Mq to Mw driven by large Yukawa couplings 
and partly due to mixing effects in the mass matrices ( see, e.g., [12] Eqns.3 - 11). As a 
result the gluino preferentially decays into third generation quarks via three body modes 
mediated by the relatively light squarks. Therefore dominant decay modes of the gluino 
are 

g-f/i + X, g^tb(tb)+Y 

where fi = t or b, X = Xi (the ith neutralino), i = 1-4 and Y= xf (the jth chargino), 
j=l,2. Hence the signals of gluino decay at the LHC will be rich in b-jets [12]. In addition 
one or more reconstructed W bosons in the final state may also help SUSY search via this 
channel. 

The size of the signal involving involving > 3 tagged b-jets (introduced mainly to 
control the QCD background) and various number of leptons was estimated by a parton 
level Monte Carlo in [12]. The SM background from tttt production only was estimated 
to be small after imposing a strong missing transverse energy Fjt cut. 

Subsequently the LHC signal of focus point supersymmetry have been studied 
by several groups using simulations beyond the naive parton level [13] - [19]. Some of 
these works required > 2 tagged b-jets in the final state for efficient background rejec- 
tion. These works will be briefly reviewed in a later section. However, none of these 
works considered the background from /1/1/2/2 production, where fi — t or b. Apriorily, 
however, some of these backgrounds with production cross sections larger than the gluino 
pair production cross section( see section 2) by several orders of magnitude seem to be 
potentially dangerous. Moreover, these backgrounds are rich in b-jets. Hence b-tagging 
which effectively suppresses light flavor QCD events may not be very effective. 

The purpose of this paper is to study the above backgrounds systematically (see 
section 2) along with a thorough study of the signal and other backgrounds. We have cal- 
culated the new background cross sections using CalcHEP[20] v2.4.5, ALPGEN [21] v2.1.1 
and MadGraph [22] v4.1. The results agree fairly well. The parton level events are then 
interfaced with Pythia [23] v6.4.8 for implementing hadronization, fragmentation and a 
toy detector simulation and analyzed with suitable selection criteria. We shall also study 
the response of these backgrounds to the cuts used by other groups [18, 19] . We next 
study the dominant QCD and tt backgrounds as well as the signal using our selection 
criteria and compare and contrast the results with other published works ( see section 2 ). 
In order to facilitate the comparison we required like other groups > 2 b-jet tags. Finally 
the consequences of selecting a signal with a richer b-jet content as suggested in [12] are 
examined. Our conclusions and future outlooks will be summarized in section 3. 
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2 Signals of Focus Point SUSY at LHC and the Back- 
grounds 



We begin by briefly reviewing some of the earlier works. Baer et al. [13] scanned the focus 
point region of the mSUGRA parameter space consistent with the WMAP data . The 
LHC reach in channels with m-jets + n-leptons + Fjt was then examined. It was concluded 
that the part of the focus point region corresponding to gluino masses < 1.8 TeV is within 
the reach of LHC. The suppression of the SM background due to light flavor QCD, tt etc 
was considered in estimating the mass reach. However, one of the key features of the 
signal in the focus point region, namely the large number of b-jets in the final state was 
not utilized in improving the signal to background ratio. 

Subsequently Mercadante et al. [15] estimated that requiring final states with 
tagged b-jets as suggested in [12] and by adjusting the kinematical cuts the gluino mass 
reach can be improved by about 20%. However, the potential backgrounds from four 
heavy flavor production mentioned in the introduction were not included in either of the 
above analyses. The viability of SUSY signals in the regions of the parameter space 
consistent with dark matter data was also studied by the CMS collaboration [16]. This 
analysis included only a part of the FP region (m < 2 TeV) and the above potential 
backgrounds were neglected. 

A few points (FP1 - FP5) in the FP region allowed by the WMAP data from 
[17] are shown in Table 1. All masses and mass parameters in this paper are in GeV 
unless otherwise stated. The other mSUGRA parameters are A = 0, tan (3 = 30 and 
sign(yu)> 0. For each point the gluino mass and the lowest order gluino pair production 
cross section (in fb) for two choices of the QCD scale (see the last two columns in Table 
1) are also shown. The first choice yielding smaller cross sections will be subsequently 
referred to as the conservative choice. The second choice, motivated by the fact that at 
this scale the leading and next to leading order production cross sections almost agree 
[24] , is the optimistic choice. It is to be noted that we have considered gluino masses near 
the LHC search limit only. The gluino pair production cross section in the leading order 
is computed by CalcHEP. The expected sharp fall of the cross section for gluino masses 
near the kinematic limit is quite apparent. 

From the dominant gluino decay modes discussed in the introduction it is clear 
that a typical final state will consist of four third generation quarks and their decay 
products. Additional quarks, leptons and missing energy would come from the decays of 
heavier charginos and neutralinos also produced in gluino decays. Thus the multiplicity 
of jets and hard isolated leptons in the signal would be large. Keeping this in mind the 
criteria listed below for selecting ( rejecting) the signal ( the background) events have 
been formulated. It should be noted that instead of separately requiring large number of 
jets and hard isolated leptons ( henceforth jets and isolated leptons will be collectively 
called objects ) we require the total number objects in the final state to exceed nine. 

• N jet > 6 and N obj > 9. 
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• Nbtag > 2 and one of the tagged jets with Et > 300 GeV. 

• Er > 300 GeV where Er has been computed using all visible objects in an event. 

• M cff = Er + Eobj E T > 2000 GeV. 

It may be noted that at this stage the constraint on Nbtag is similar to the one 
employed by other analyses. We shall show below that requiring more tagged b-jets in 
the final state as suggested in [12] improves the signal vis a vis the background. 

We begin the analysis by estimating the backgrounds from production of four 
heavy flavors in different combinations. This will be followed by the simulation of more 
conventional backgrounds and the signal also studied by other groups. The production 
cross sections computed by CalcHEP, ALPGEN and MadGraph are presented in Table 2. We 
have set the following input values at the weak scale (M wea k): a cm = 127 1 934 , a s =0.1172, 
m t (pole)=175.0 [25], m b (m b ) = 4.25 GeV. 

In our calculation we use the CTEQ5L [27] parton distribution functions (PDF). 
The appropriate loop corrected values of strong coupling constant at the renormalization 
scale is calculated on the fly during the convolution with the PDF. To avoid the generation 
of unwanted events the following nominal cuts are applied only on the bottom quarks: 
p T >10 GeV, 77(b) <5.0 and AR(b,b) >0.3 where AR = ^A?] 2 + A0 2 . 



QCD scale 



FP point ( m , m 1/2 ) 


m i 




m i/2 


FP1 ( 3700.0, 700.0 ) 


1751.3 


1.18 


2.37 


FP2 ( 3975.0, 730.0 ) 


1824.7 


0.81 


1.63 


FP3 ( 3975.0, 790.0 ) 


1950.0 


0.39 


0.80 


FP4 ( 4130.0, 833.0 ) 


2047.4 


0.24 


0.48 


FP5 ( 4235.0, 866.0 ) 


2119.8 


0.16 


0.33 



Table 1: Cross-sections (in fb) for pp — > gg at LHC using CalcHEP v2.4.5 for different 
FP scenarios and two choices of the QCD scale (see the last two columns). The masses 
and mass parameters are in GeV. 

In Table 2 we demonstrate the variation of the cross sections with the QCD scale. 
There is reasonable agreement between different generators at the same scale. As expected 
there is significant scale dependence in these leading order cross sections. However, in 
spite of this uncertainty some of the background processes have cross sections several 
orders of magnitude larger than that of the signal, which is typically a few fb for m~ near 
the LHC search limit (see Table 1). Thus it is worthwhile to see the response of these 
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Process Generator 0.5^ 2.0^ 



tttt MadGraph 4.38 2.82 1.93 

ALPGEN 4.42 2.90 1.96 

CalcHEP - 2.89 

ttbb MadGraph 2.62 x 10 3 1.76 x 10 3 1.14 x 10 3 

ALPGEN 2.86 x 10 3 1.89 x 10 3 1.29 x 10 3 

CalcHEP - 2.27 x 10 3 

bbbb MadGraph 1.20 x 10 7 7.61 x 10 6 5.62 x 10 6 

ALPGEN 1.22 x 10 7 8.04 x 10 6 5.52 x 10 6 

CalcHEP - 1.15 x 10 7 



Table 2: Production cross-section (fb) for different final states involving four third gen- 
eration quarks obtained using different generators. Variation of the cross sections with 
QCD scale is shown in the last three columns. 

backgrounds to our choice of the kinematical cuts listed below as well as to the cuts used 
by similar analyses [18, 19]. We shall subsequently show that our cuts also lead to a 
healthy signal size and control other conventional backgrounds quite effectively. 

We next simulate the backgrounds listed in Table 2. The parton level unweighted 
events from ALPGEN has been interfaced with Pythia. As noted earlier our final conclusion 
regarding the importance of these backgrounds remain unchanged inspite of the sizable 
uncertainty due to the choice of scale. In order to avoid generation of events unlikely to 
pass the cuts we have imposed the following conditions on the parton level sub-processes. 
It should be noted that these conditions are somewhat stronger than those used in Table 
2. 

• tttt events have been generated with no cut on t. 

• ttbb events have been generated with no cut on t, P^ > 20 GeV and \rjb\ < 4.5. 

• bbbb events have been generated with P^ > 50 GeV and \r] h \ < 4.5. 

We have used the toy calorimeter simulation (PYCELL) provided in Pythia with 
the following criteria: 

• The calorimeter coverage is < 4.5. The segmentation is given by At] x A0 = 
0.09 x 0.09 which resembles a generic LHC detector. 
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• A cone algorithm with AR = yj At] 2 + A<p' 2 = 0.5 has been used for jet finding. 

• E J t min = 30 GeV and jets are ordered in E T . 

• Leptons ( I = e, /i ) are selected with Et > 20 GeV and \q\ < 2.5. 

• No jet should match with a hard lepton in the event. 

• A jet with \r}\ < 2.5 matching with a B-hadron of significant decay length has been 
marked tagged. 

We have implemented jet and lepton( i = e or /x) isolation using the following 
criteria: if there is a jet within AR = 0.5 and E^/E^ < 1.2, the jet is removed from 
the list of jets, else the lepton is removed from the list of leptons. Single b-jet tagging 
efficiency in tt has been tuned to eb — 0.5. 




Missing E T (GeV) K« (GeV) 



Figure 1: Distributions for missing transverse energy (Ep) (left) and effective mass (M e ff) 
(right) for four heavy flavor final states and signal with m~ = 1751.3 GeV ( FP1, see 
text): in both the figures the distributions are indicative of shapes only and normalized 
to unit cross section. 

The distributions of Er and M e g for different final states involving four third gen- 
eration quarks are shown in Fig. 1. Selection efficiencies for these background processes 
due to individual cuts are shown in Table 3; the individual cut which rejects most events 
for each final state is shown in bold. The residual cross-sections for these backgrounds 
after all our selection cuts are shown in Table 4. As is expected from Fig. 1 very little 
background survives. 

For comparison we have computed the backgrounds corresponding to the cuts 
implemented by the ATLAS collaboration [18] (column 5 of Table 4) and Baer et al. [19] 
(column 6 of Table 4) . It may be noted that in all cases these backgrounds are rather 
small compared to other backgrounds from tt production or QCD processes (see below). 
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Process N jct > 6 N obj > 9 N btag > 2 Pf^ 1 > 300 Br > 300 M eff > 

GeV GeV 2000 GeV 



tttt_ 0.9304 0.5457 0.5946 0.1014 0.0316 0.0346 
ttbb_ 0.5792 0.0831 0.4509 0.0265 0.0033 0.0018 

bbbb 0.2186 0.0064 0.5854 0.0120 0.00004 0.0002 



Table 3: Selection efficiencies for different final states involving four third generation 
quarks due to individual cuts have been shown. The individual cut which rejects most 
events for each background is shown in bold. 



Process 


Cgcn 


Events 


This work 


ATLAS 


Baer et al. 




(fb) 


generated 


^sel • ^gcn 


^sel • Cgen 


^sel • "gen 


tttt 


5.3 


1.25 x 10 4 


0.012 


0.093 


0.004 


ttbb 


831 


1.40 x 10 5 


0.042 


0.665 


0.059 


bbbb 


7260 


3.02 x 10 5 


0.000 


0.240 


0.000 



Table 4: The cross section at the generation level (o" gcn ), number of events generated, and 
the residual cross-section (e sc i.o" gen ) after all cuts. The cross sections were computed with 
cuts on final states partons different from the ones used in Table 2 and the default QCD 
scale ( Q 2 = S f (m 2 + P^), where '¥ stands for the final states partons) of ALPGEN. 
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This conclusion holds even if a factor of two due to the QCD scale uncertainty is taken into 
account. It may be recalled that the cuts in [18] are designed for relatively light gluinos in 
the focus point region. Consequently they are less severe than the ones introduced in this 
paper or in [19]. As a result some of the residual backgrounds of this type are numerically 
significant for integrated luminosities > 100 fb _1 and require some attention. 

Next we turn our attention to the conventional backgrounds. First we shall 
summarize the earlier works which employed b-tagging. Mercadante et al. [15] found that 
the dominant background comes from QCD processes involving light flavors and gluons 
only. The analysis of reference [18] ignored this background arguing qualitatively that 
the requirement of 2-tagged b-jets should be adequate to reduce it to a negligible level. 
However, they studied the bbj background, which is included in the QCD background in 
this work as well as in [19] and found it to be 41% of the dominant tt background. Very 
recently Baer et al. [19] analyzed the same background with different kinematical cuts. 
They found that this background is rather small. 

It has been identified that the tt background is the dominant one in both [19] 
and [18] and we agree with them qualitatively. In the former work this background was 
simulated by ISA JET [28]. However, in addition to other selection criteria they required 
> 7 jets + > 2 tagged b-jets in the final state. Since the parton level process has at 
most six jets (including two b-jets) the size of this background crucially depends on the 
model of parton showering in the generator. Our cut on N Q bj is similar in nature. We 
shall briefly comment on the reliability of analyses requiring high jet multiplicity in the 
following. 

We now present our estimates of these backgrounds. The cross-sections for QCD 
processes (which includes bb and cc production) and tt production are quite large (see 
Table 6). Since the selection demands a large number of jets, Br and M eff it is more likely 
that events with higher will survive. Even then both the above backgrounds survive 
the cuts (mainly the one on the multiplicity of objects) because of initial state and final 
state radiation (ISR and FSR) and multiple interaction. On the other hand events with 
lower y/s have a very large share of the cross-section and ISR and FSR help them to pass 
the selection. To estimate the contributions from the above two regions with sufficient 
statistics and to avoid generating large number of events which will not pass the selection, 
we have generated events with certain cuts at the parton level: 

• tt events have been generated in two bins, < 600 GeV (tt-low) and P^ > 600 
GeV (tt-high). 

• QCD events have also been generated in two bins, 100 GeV < Pt < 600 GeV 
(QCD-low) and P T > 600 GeV ( QCD-high) for the final state partons. 

The efficiency of the individual cuts and the size of the different backgrounds are 
listed separately in Tables 5 and 6. The cross sections are calculated with the conservative 
choice of the QCD scale. It is to be noted that while the cut N obj > 9 is most effective for 
rejecting tt-high events with high y/s, it is either the cut on Br or M eff for events with low 
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y/s. The Fir distributions before and after applying all other kinematical cuts are shown 
in Fig. 2. The distributions of Nj et and N Q bj are shown in Fig. 3. It is obvious that 
even after the cut Nj et > 7 quite a bit of background survives where as the cut N Q bj > 9 
suppresses the background more effectively. 




200 400 600 800 1000 200 400 600 800 1000 

Missing E T (GeV) Missing E T (GeV) 



Figure 2: Missing transverse energy (Fjt) distributions for the major backgrounds and the 
signal (m~ = 1751.3 GeV). In the left plot the distributions are indicative of the shapes 
only. In the right plot they have been normalized to the residual cross-sections after all 
cuts except the one on Fir- 

We do not impose any extra strong hardness cut on a jet which is counted as an 
acceptable object. It is quite plausible that in a real life situation such jets indeed arise 
from showering. However as noted before the number of such jets in an event crucially 
depends on the parton showering model in the generator. Therefore the reliability of the 
background estimate both in this work as well as in [19] can not be guaranteed until the 
showering model is validated using LHC data. In the next paragraph we briefly summarize 
a study [29] which indicates that the tt background estimated as above may not be wide 
off the mark. 

It is certainly desirable to supplement the lowest order estimates of the back- 
ground by higher order matrix element calculations. However, care must be taken to 
avoid double counting of jets which may arise in a lower order calculation due to hard 
emissions during shower evolution as well as from a higher order matrix element. A proper 
matching of matrix elements and shower evolution is, therefore, called for. In reference [29] 
the lowest order cross section of tt production was computed. Next the combined cross 
section of the processes tt (exclusive), tij (exclusive), tijj (exclusive) and tijjj (inclusive) , 
where j stands for a light quark or gluon, were obtained after matrix element matching( 
see Table 1 and Table 2 of [29]). The two results agree nicely both for Tevatron and LHC 
energies. A similar agreement is found for various distributions of the produced tt system 
computed by the two different ways. It is these distributions which after all determine 
the efficiencies of various cuts. Thus one may conclude that at least the lowest order tt 
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Figure 3: The distributions of the number of jets (Nj et ) (top) and number of objects (N Q bj) 
(bottom) for the major backgrounds and the signal (m~ = 1751.3 GeV). The figures on 
the left show the distributions before any cut (normalized to unit cross section); the figures 
on the right show distributions normalized to residual cross-sections after all cuts except 
those on Nj et (top) and N ob j (bottom). 
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Process 


Niet > 6 


Nobi > 9 


Nbtae > 2 


P^ 1 > 300 


Br > 300 


M eff > 










GeV 


GeV 


2000 GeV 


Signal 


0.9674 


0.7116 


0.7899 


0.6275 


0.7416 


0.8173 


QCD - high 


0.1490 


0.0074 


0.0177 


0.0599 


0.0018 


0.1033 


QCD - low 


0.0151 


0.0001 


0.0052 


0.0003 


0.0000 


1.9 x 10~ 6 


tt - high 


0.2318 


0.0188 


0.4567 


0.5853 


0.1799 


0.0963 


tt - low 


0.1984 


0.0083 


0.1490 


0.0074 


0.0023 


0.0014 



Table 5: Selection efficiencies for major backgrounds due to individual cuts. The individ- 
ual cut which rejects most events for each case is shown in bold. 



Process 


CTg Cn 


Events Residual cross-section (e sc i.o- g cn) 




(fb) 


generated N btag > 2 N btag > 3 



Signal 

m i = 1751.3 
mj = 1824.7 
= 1950.0 
= 2047.4 
= 2119.8 



m- 
m- 



1.20 
0.81 
0.40 
0.24 
0.16 



10 4 
10 4 
10 4 
10 4 
10 4 



0.3598 
0.2267 
0.1445 
0.0921 
0.0652 



0.2445 
0.1527 
0.0987 
0.0634 
0.0431 



Background 
QCD - high 
tt - high 
tt - low 



1.81 x 10 5 
4.14 x 10 2 
3.70 x 10 5 



1.08 x 10 7 
2.325 x 10 6 
1.14 x 10 7 



0.4362 
0.1675 
0.8806 



0.0504 
0.0205 
0.1793 



Table 6: The cross-section at the generation level (o" gcn ) , number of events generated, 
and residual cross-section (e se i.cr ge n) after all cuts are shown. <r gcn is calculated for Q 2 = s 
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signal 
QCD - high 
tt - high 
tt — low 



0.7 
0.6 
0.5 



b 
-I 



b 
- I 



0.3 
0.2 
0.1 




Figure 4: Distributions for number of b-tagged jets ( N btag ) are shown: (left) before 
any selection criteria has been applied and (right) after other selection criteria have been 
applied. In both the distributions are indicative of shapes only. 

background including parton showering which is the dominant background for the signal 
under study is fairly reliable. 

The case of the pure light flavor QCD background is rather similar. The required 
final state with nine objects from this 2^2 process can arise solely due to parton 
showering and, therefore, model dependent to some extent. Simulating the lowest order 
QCD processes in Pythia we find that the size of this background is smaller than the tt 
background but is non-negligible. However the entire contribution comes from QCD-high 
events P T > 600 GeV. After generating 15 million QCD-low (100 GeV < P T < 600 GeV) 
events we found that no event survive the cuts. 

We now compare our results with other calculations which required two or more 
tagged b-jets in the final state. A comparison with [18] may not be meaningful as their 
cuts are designed for the search of a relatively light gluino (m~ = 857.0). In fact with 
their cuts we find that the survival probability of the tt background is twice as large as 
our estimate and will overwhelm the gluino signal for m~ pa 2 TeV, which is roughly the 
LHC reach in the FP region, having a cross section not exceeding a few fb's. 

The tt background after all cuts estimated by Baer et al. [19] is 1.5 fb with 'a tiny' 
contribution from QCD. On the other hand we get a smaller tt background while the QCD 
background is non- negligible. Thus one may conclude that our cuts are more efficient in 
reducing the former background while the cuts of [18] may be better for eliminating the 
QCD backgrounds. A more detailed comparison is not possible since we donot know the 
QCD scales for different cross sections in [19]. 

It is probably a numerical conspiracy that the total background (tt + QCD) 
computed by the two groups agree rather nicely. However, the fact that two groups using 
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different event generators and different sets of kinematical cuts get comparable results at 
least indicates that the major backgrounds obscuring the search for focus point SUSY at 
LHC are fairly well understood. 

We next present the signal size for different gluino masses. The efficiencies of the 
cuts used in this paper for the scenario FP1 are in Table 5 while the signal size for the 
conservative scale choice in different scenarios after all cuts are in Table 6. The signal 
cross sections should be multiplied by a factor of 2 if the scale m~/2 is considered (see 
Table 1). The Bp, Nj et and N Q bj distributions are in Figs. 2 and 3. Comparing these with 
the background cross sections and taking the QCD scale uncertainty into account it is 
clear that the gluino mass reach at the LHC cannot be pin pointed as yet. Nevertheless 
the following comments are in order. 

Using the optimistic scale choice (= m~/2) one finds that for FP1(FP3) (corre- 
sponding to m~ = 1751.3 (1950)) S/V^IO.O (4.0) with an integrated luminosity 300 
fb _1 , where S and B are the number of signal and background events respectively. Of 
course the S/v^B ratio may be modestly increased by optimizing the cuts [15]. On the 
other hand the LO signal cross section at the scale m~/2 mimics the NLO cross section 
[24]. In view of this the next to leading order ti cross section (about a factor of 2 larger 
than the cross section in Table 6 ) may give a more reliable background estimate. This 
will further suppress the S/\^B ratio by a factor of \f2. For the conservative scale choice 
on the other hand the above ratios will be suppressed by a factor of two. On the whole 
it seems unlikely that using the selection criteria based on 2 or more tagged b-jets as is 
usually required, the gluino mass reach at LHC can be pushed far beyond 1.8 TeV with 
certainty. 

We now impose a more stringent b-identification criteria as suggested in [12] and 
require > 3 b-jet tags keeping all other selection criteria listed above unchanged. The 
residual signal and the backgrounds are shown in the last column of Table 6. Even the 
conservative scale choice now yields S/v 7 ^ ratio 8.5 (3.4) for m^ = 1751.3 (1950). For the 
optimistic scale choice the above ratios are enhanced approximately by a factor of V2.0 
after doubling the background which should reflect the theoretical uncertainties due to 
possible NLO corrections etc. In view of the possibility of further improving this ratio by 
optimizing the selection criteria, one can hope that irrespective of the uncertainty due to 
the choice of scale the reach in m~ at LHC is close to 2.0 TeV. 



3 Conclusion 

If supersymmetry in the focus point region is realized in nature it may lead to spectacular 
signals at the LHC via gluino pair production [12, 13, 15, 16, 17, 18, 19]. We propose a 
new set of criteria for selecting(rejecting) the signal( background) (see section 2). 

We have extended the estimation of the backgrounds which may hinder the search 
for supersymmetry in the focus point region. In particular we have simulated and ana- 
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lyzed the potentially dangerous backgrounds from different processes involving four heavy 
flavor quark production which have so far been neglected. Some of these processes have 
cross sections several orders of magnitude larger than that of the signal. Moreover these 
backgrounds are rich in b-jets and cannot be reduced by b-tagging which is recommended 
as one of the most effective tools for removing other SM backgrounds [12, 15, 19]. 

Using rather stringent cuts (see section 2) appropriate for detecting relatively 
heavy gluinos with masses near the kinematic reach of the LHC, we have shown that these 
background are highly suppressed(see Table 4). On the other hand for less severe cuts 
often employed for relatively light gluinos in the focus point region (see for example [18]), 
these backgrounds may be small but numerically significant and call for some attention. 

We have also simulated the conventional SM backgrounds estimated by other 
groups. We qualitatively agree with the observation [18, 19] that tt events are the most 
severe background to the gluino signal. In view of various theoretical uncertainties in 
the production cross sections the LHC reach in m~ cannot be pinpointed as yet. This 
is especially so if the selection criteria is restricted to > 2 tagged b-jets which is the 
usual strategy. Requiring 3 or more tagged b-jets in the signal the gluino mass reach at 
LHC is close to 2.0 TeV even with a modest production cross section corresponding to a 
conservative choice of the QCD scale. 

We look forward to improved b-tagging efficiencies along with improvement in 
the theoretical understanding (e.g., the inclusion of matrix element matching) of the light 
flavor QCD and other backgrounds for a better handle on gluino searches in the FP 
scenario. 
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